Excluded Volume Effects in the Quark Meson Coupling Model 
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Excluded volume effects are incorporated in the quark meson coupling model to take into account 
in a phenomenological way the hard core repulsion of the nuclear force. The formalism employed is 
thermodynamically consistent and does not violate causality. The effects of the excluded volume on 
in-medium nucleon properties and the nuclear matter equation of state are investigated as a function 
of the size of the hard core. It is found that in-medium nucleon properties are not altered significantly 
by the excluded volume, even for large hard core radii, and the equation of state becomes stiffer as 
the size of the hard core increases. 



I. INTRODUCTION 



The study of the properties of high density and high temperature hadronic matter is of interest for understanding a 
wide range of phenomena associated with superdense stars |Q and relativistic heavy-ion collisions . One of the open 
questions in this subject is the correct identification of the appropriate degrees of freedom to describe the different 
phases of hadronic matter. Although this question will eventually be answered with first-principles calculations with 
the fundamental theory of the strong interactions quantum chromodynamics (QCD), most probably through lattice 
QCD simulations, presently one is still far from this goal and, in order to make progress, one must rely on model 
calculations and make use of the scarce experimental information available. For matter at zero temperature and 
density close to the saturation density of nuclear matter, experiment seems to indicate that the relevant degrees 
^ of freedom are the baryons and mesons. There is a long and successful history of calculations using models based 
on baryonic and mesonic degrees of freedom, such as potential models ||[^] and relativistic field-theoretical models, 
generically known as quantum hadrodynamics (QHD) [^|J^]. At densities several times larger than the saturation 
density and/or high temperatures, one expects a phase of deconfmed matter whose properties are determined by the 
li ■ internal degrees of freedom of the hadrons. Early studies of deconfmed matter used the MIT bag model Jg], in 
which the relevant degrees of freedom are quarks and gluons confined by the vacuum pressure. On the other hand, at 
high densities, but not asymptotically higher than the saturation density, the situation seems to be very complicated. 
The complication arises because of the possibility of simultaneous presence of hadrons and deconfmed quarks and 
gluons in the system. Not much progress has been possible in this direction because of the necessity of a model able 
to describe composites and constituents at the same footing - a step towards this direction is the formalism developed 
in Ref. §. 

An important step towards the formulation of a model to describe the different phases of hadronic matter in terms 
of explicit quark-gluon degrees of freedom is the quark-meson coupling (QMC) model, proposed by Guichon pi 
some time ago and extensively investigated by Saito, Thomas and collaborators - see also Ref. Q for related 
studies. Matter at low density and temperature is described as a system of nonoverlapping bags interacting through 
effective scalar- and vector- meson degrees of freedom, very much in the same way as in QHD (^||. The crucial 
difference is that in the QMC, the effective mesons couple directly to the quarks in the interior of the baryons, with 
the consequence that the effective baryon-meson coupling constants become density dependent. In addition, hadronic 
sizes are explicitly incorporated through form factors calculable within the underlying quark model Jl3[ ] . At very high 
density and/or temperature, baryons and mesons dissolve and the entire system of deconfmed matter, composed by 
quarks and gluons, becomes confined within a single MIT bag @. 

The fact that the same underlying quark model is used at different phases of hadronic matter makes the model 
very attractive conceptually. Many applications and extensions of the model have been made in the last years - see 
Refs. |l3|-|l9| and references therein. Of particular interest for the phenomenology of finite nuclei was the introduction 
of a density-dependent bag constant by Jin and Jennings |p^| . These authors postulated different density dependencies 
for the bag constant, in a way that the bag constant decreases as the nuclear density increases. One consequence 
of this is that large values for the scalar and vector mean fields at the saturation density are obtained, leading to 
spin-orbit splittings of the single-particle levels of finite nuclei that are in better agreement with experiment than 
those obtained with a density independent bag constant. Another consequence of a smaller bag constant in medium 
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is that the bag radius becomes considerably larger than in free space Jl4|,[l6| . At saturation density the nucleon radius 
increases 25% and at densities three times higher the radius can increase as much as 50%. For higher densities the 
increase of radius is even more dramatic. The consequences of changing the bag constant for nucleon sizes were 
investigated by Lu et al. 

A large increase of the bag radius naturally raises the question about the validity of the nonoverlapping bag picture 
that underlies calculations of nuclear matter properties with the model. At normal nuclear matter densities, the 
average distance between nucleons is of the order of 1.8 fm. Therefore, for densities larger than the normal density and 
bag radii larger than 1 fm there is a large probability that the bags overlap significantly. However, before concluding 
that the picture of independent nucleons breaks down, it is important to recall that short-distance correlations are 
left out in a mean field calculation. These correlations are induced by the combined effects of the Pauli exclusion 
principle between identical nucleons and the hard-core of the nucleon-nucleon interaction that forbids scattering into 
occupied levels. The success of the independent particle model of the nuclear shell model is due to the small size of 
hard core and the Pauli Principle that lead to a "healing" distance of the two nucleon relative wave function that 
is smaller than the average distance between nucleons in medium [^(J. In a model like the QMC, where the finite 
size of the nucleons is made explicit through a bag structure, the incorporation of this physics in the many-body 
dynamics is an interesting new development. In the present paper we address this in a phenomenological way through 
an excluded volume approach. The prescription we use was developed in Ref. [^lf for ideal gases and further extended 
to relativistic field-theoretic models as QHD in Ref. p2| . In this approach, matter in the hadronic phase is described 
by nonoverlapping rigid spheres, but when the density of matter is such that the relative distance between two spheres 
becomes smaller than the diameter of the spheres, the excluded volume effect introduces an effective repulsion that 
mimics the hard-core repulsion of the nucleon-nucleon interaction. 

Of course, at very high density the description of hadronic matter in terms of nonoverlapping bags should break 
down. In a purely geometrical view, one has the picture that once the relative distance between two bags becomes 
much smaller than the diameter of a bag, quarks and gluons start percolate and individual bags loose their identity. 
The density for which this starts to happen is presently unknown within QCD. In this respect, it is important not 



to confuse the bag radius with the radius of the hard-core of the nucleon-nucleon force. Model studies 23 indicate 
that when the two nucleons start to overlap, medium-range forces are generated from the distortion of the quark 
distribution. The short-range repulsion, on the other hand, is due to the combined effect of the one-gluon exchange - 
mainly due to its spin-spin component - and the Pauli exclusion principle between quarks of different nucleons that 
becomes efficient when the overlap of the two-nucleon wave functions is complete. Although the described scenario 
might well not be ultimately confirmed by a full QCD calculation, it seems however clear that the two radii, the 
radius of the MIT bag and the radius of the hard-core of the nucleon-nucleon interaction, are of different sizes and 
have different origin in the physics of hadron structure. In this sense, the volume of the excluded volume will be taken 
to be smaller and unrelated to the radius of the underlying MIT bag. 

The excluded volume approach we use is thermodynamically consistent. Although the prescription can be extended 
to take into account Lorentz contraction of the bags |24|] , in this initial exploratory investigation we use hard-sphere 
bags, since a complete calculation would lead to massive numerical calculations. In addition, as indicated by the 
investigations in Ref. pi| , the effect of Lorentz contraction is most important for light particles like pions. However, 
as we will explicitly show, the approach does not lead to violations of causality for the density range where the model 
makes sense. 

The paper is organized as follows. In Section || we present a short review of the excluded volume prescription of 



Ref. |22J and implement it to the QMC model. Numerical results are presented in Section III and our Conclusions 



and Perspectives are discussed in Section [V 



II. EXCLUDED VOLUME IN THE QMC MODEL 

Initially, for completeness and in order to make the paper selfcontained, we briefly recapitulate the excluded volume 
prescription of Ref. [ p2[ . Let us start with the ideal gas of one particle species with temperature T, chemical potential 
/j, and volume V. The pressure is related to the grand partition function Z as 

P(T,,)= lu*T lnZ ^ V \ (1) 

with Z defined as 

oc 

Z(T, H,V) = J2 ^ N/T Z{T, N, V), (2) 

N=0 
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where Z is the canonical partition function. The authors in Ref. |2^] included the excluded volume effect starting 
from the canonical partition function as 



zexcl 



(T, N, V) = Z(T, N,V — v Q N)S(V - v N). (3) 



This ansatz is motivated by considering the volume V for a system of N particles is reduced to an effective volume, 
V — vqN where vq is the volume of a particle. In an hadronic gas, vq can be interpreted as the region excluded by the 
hard core of the nuclcon-nucleon interaction. For a spherical region, vo — 47rr 3 /3 with r the hard-core radius. Using 
Eq. (|J) into Eq. (||) , the grand partition function becomes 

oo 

Z excl (T, fj,,V) = J2 e-» N/T Z(T, N,V — v N)e(V - v N). (4) 

There is a difhculty for evaluation of the sum over N particles in this equation because of the dependence of the 
available volume on the varying number of particles N because Z(T, N, V — VqN) does not factor as a product as in 
the case of an iV-independent volume. To overcome this difficulty the authors in Ref. 22 have performed a Laplace 
transformation on the variable V in Eq. (^) as 

POO 

Z excl (T, dVe~^ v Z excl {T,p,V). (5) 

Jo 

Using Eq. (^) in this, and making the change of variable 

V = x + v Q N, (6) 

one obtains 

/•oo 

Z™ cl (T,n,0= dxe-^Z(T,fi,x)=Z excl (T,fi,0, (7) 
Jo 

where p, = fj, — VqT^. Now the integrand in this is factorizable (for the present case of independent particles) and the 



sum over N can be implemented. It is a simple exercise ElJ to show that the pressure of the system is given as 22 



P(T, M ) = P'(T,/i), (8) 

with 



fl = fx-v P(T,ix). (9) 

The meaning of Eq. (||) is that the pressure of the system with excluded volume with chemical potential /i, P(T, fi), is 
equal to the pressure of a system without excluded volume but with an effective chemical potential ft = \i — v$P(T, /x), 
denoted by P'(T,p). Note that once the expression for P'(T,ji) is known, the pressure of the system is given by an 
implicit function. 

The baryon density, the entropy density and the energy density for the system are given by the usual thermody- 
namical expressions 

P(T, M) 
5(T )M ) 

These relations define a thermodynamical consistent formalism, since the fundamental thermodynamical relations are 
fulfilled. 

Next we apply this formalism to the QMC model for nuclear matter at zero temperature In the QMC 

model, the nucleoli in nuclear matter is assumed to be described by a static MIT bag in which quarks interact with 
scalar a and vector cu mean mesonic fields. The mesonic fields are meant to represent effective degrees of freedom, 



dP 

>, 
' dP S 
df 



TS-P + f ip 



p'(T,ji) 
l + v (/(T, 

S'(T,~P) 
1 + v oP'(T, ft)' 



1 + V p'(T,Ll) ' 



(10) 

(11) 
(12) 
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not necessarily identified with real mesons. Therefore, since the introduction of the excluded volume is to represent 
the hard core nucleon-nucleon interaction, Eqs. ([§|)-(|l4) will be applied to the baryons only. The same prescription 
has been used in the application of the formalism to QHD in Ref. [^2| . 

In the QMC model, the pressure and energy density receive contributions from baryons and mesons and are given 

as 



P = P R - -mloZ 



1 



(13) 
(14) 



where Pb and 6 b are the baryon contributions. As said above, the excluded volume will be applied only to the 
baryonic contributions and, at zero temperature, one needs to consider only Eqs. (JTo|) and (p"2|), 



P = 



P 



1 + v oP ' 
j 



1 + v p' 

For practical calculations, it is convenient to parameterize p' in terms of a kp according to 

o> - -L k 3 

This allows to write the QMC expressions for P' B and e' B as 



(15) 
(16) 

(17) 



P'b 



3 2vr 2 



kj?\/ k 2 w 



M* 



M* z k F \ kl + M* 



+ ?M* 4 ln 



kp 



k\ + M* 



M* 



e' B = P J 



kp 



M 



*2 



P'b- 



(18) 



(19) 



In these, M* is the in medium nucleon mass calculated with the MIT bag. Its value is determined by solving the 
MIT bag equations for quarks coupled to the mean fields gq and ujq. In order to completely determine M* , and 
therefore the nuclear matter properties Eqs. ( |l5| ) and (|l6|), one needs <jq and ujq. The scalar mean field is determined 
selfconsistently from the minimization condition at density p: 



de 
da 



= 0, 



which leads to 



with 



1 



EM 



1 + v p' ml 



1 dM* 



k F \ kl + M* 2 -M* 2 ln 



' k F + ^k 2 F + M* 
M* 



The vector mean field u>o is obtained from its equation of motion as 



OJQ 



■P- 



(20) 



(21) 



(22) 



(23) 



Solution of Eq. ( pl| ) proceeds as follows. For a given p, we use Eq. (|T^) to obtain p' , and form this p' we obtain kp 
of Eq. (|l7|). The derivative dM* /duo can be done explicitly, S(cto) is then known, and the transcendental equation 
for Co is easily solved numerically The results are presented in the next section. 



4 



III. RESULTS AND DISCUSSIONS 



We start fixing the free-space bag properties. We use zero quark masses only and use two values for the bag radius, 
R = 0.6 fm and R = 0.8 fm. There are two unknowns, zq and the bag constant B. These are obtained as usual by 
fitting the nucleon mass M — 939 MeV and enforcing the stability condition for the bag. The values obtained for zq 
and B are displayed in Table 1. 

Next we proceed to nuclear matter properties. We will consider two versions of the model. In the first one, the bag 
constant B is fixed at its vacuum value, and in the second one the bag constant changes accordingly to the original 
Jin and Jennings |3 ansatz, namely 

B* = Bc X p(-^-), (24) 



N 



where g^ is an additional parameter and B is the value of the bag constant in vacuum. In this work we use g^ = 2.8, 
which is the same as in Ref. |l7j . 

The quark-meson coupling constants g q and <? w = 3<?2, are fitted to obtain the correct saturation properties of 
nuclear matter, Eb = E/A — M = e/p — M = —15.7 MeV at p = po = 0.15 fm~ 3 . We take the standard values for 
the meson masses, m a = 550 MeV and to w = 783 MeV. We present results for three different values of the hard core, 
r = 0.4 fm, 0.5 fm and 0.6 fm, and for two values of bag radii, R — 0.6 fm and 0.8 fm. The pair r = 0.6 fm, R = 0.6 fm 
represents the situation that the size of the hard-core is the same as of the bag and is included for illustrative purposes. 

Initially, we investigate the effect of the excluded volume on the binding energy per particle for the values of r 
and R mentioned above. The results for the different values of r and R and shown in Fig. [I], where we plot Eb as 
a function of the nuclear density p. In this figure the coupling constants g% and g u for a given R are the same for 
the different values of r. As expected, the effect of an effective repulsion due to the hard core is clearly seen in this 
figure. The effect obviously increases as the size of hard-core increases. At the saturation density, the largest value 
of the effective repulsion is of the order of 4 MeV. The effect is not as dramatic as one could expect. For comparison 
with another repulsive effect, we mention that Fock terms |l3j give 5 MeV repulsion for the binding energy. 

We now readjust the coupling constants g q a and g^ such as to obtain the correct saturation binding energy of 
nuclear matter for the different values of r and R. Our aim is to investigate the changes on the properties of nuclear 
matter and in-medium nucleon properties due to the hard core. Tables II and III present the values of the coupling 
constants and the ratios of in-medium to free-space bag radii R* /R, nucleon masses M* jM and bag eigenvalues x* fx. 
The Tables also show the changes in the incompressibility for different hard core radii. The results are such that 
nucleon properties are not changed significantly, being at most at the level of 2%. The incompressibility is a little 
more sensitive than nucleon properties to the extra repulsion induced by the hard core, but the increase is at most 
120 MeV. 

The effect of the hard core as a function of the nuclear density p on the binding energy is shown is Figure ||. One 
notices that the equation of state becomes stiffer as the size of the hard core increases. The ratios R* /R, M* jM and 
the <tq field as functions of p are shown in Figures |^, [s] and respectively. As found previously, the in-medium bag 
radius decreases (increases) for a constant (in-medium changed) bag parameter. Now, the change in the in-medium 
bag radius decreases as the hard core radius increases. This is clearly an effect due to the fact that as the hard 
core radius increases, one has less attraction, and the bag properties change less. In Fig. || one sees the interesting 
feature that as the in-medium nucleon mass increases, the binding energy curve is stiffer when volume corrections are 
included, contrary to the case without excluded volume. This is again an effect of extra repulsion due to the hard 
core. That one gets less attraction as the hard core radius increases can be seen in Figure |^, where we plot ctq as 
function of p for different combinations of r and R. The less attraction is simply due to the factor 1 + vop in the 
denominator in Eq. (pi]), which increases as vq increases and makes the r.h.s. of Eq. (|2l| ) to contribute less to do- 

To conclude this section, we mention that for neutron stars, for instance, one is interested in the equation of state 
pressure P versus energy e. One important question here is to check whether causality is respected by such an equation 
of state. Figure || presents P versus e for different values of r and R. For comparison, the causal limit P = e is also 
shown in the figure. Clearly seen is that all the cases studied here respect the causal condition dP/de < 1, so that 
the speed of sound remains lower than the speed of light. This result is consistent with Ref. p5j ], where it was shown 
that for realistic situations of temperatures below the QCD phase transition, which is believed to be of the order of 
200 MeV, the excluded volume prescription used here |2]]J22|] does not lead to conflicts with causality. 
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IV. CONCLUSIONS AND PERSPECTIVES 



In this paper we have incorporated excluded volume effects in the quark meson coupling model in a thermody- 
namically consistent manner. The excluded volume simulates in a phenomcnological way the short range hard-core 
repulsion of the nucleon-nucleon force, in the sense it does not allow nucleons to occupy all space as they were point- 
like. The consequences for in-medium nucleon properties and saturation properties of nuclear matter due to the 
excluded volume effects have been investigated for different bag and hard core radii. The bag constant was allowed 
to change in medium and differences with respect to a fixed bag constant were studied. It was also shown that the 
prescription used does not lead to violations of causality. 

We found that the excluded volume induces an effective repulsion that increases as the size of hard-core increases. 
The repulsion is at most 4 MeV at the saturation density. In- medium nucleon properties, such as bag radius and 
nucleon mass are not changed significantly, as compared to the changes when excluded volume effects are not taken 
into account. The changes are at most at the level of 2%. The incompressibility is a little more sensitive, but the 
increase is at most 120 MeV. The excluded volume also induces the effect that the binding energy curve as a function 
of the nuclear density is stiffer as the in-medium nucleon mass increases. This feature is contrary to the case without 
excluded volume. It arises because of the extra repulsion due to the hard core that leads to a smaller sigma field and 
consequently to less attraction. 

The formalism of the present paper can be extended in several ways. We intend to incorporate taking into account 
of the Lorentz contraction of the bags. As indicated by the authors of Ref. ]24[] , the effects of the Lorentz contraction 
is most important for light particles like pions. Also possible extensions of the formalism presented here to finite 
temperatures are currently under progress. As a final remark, it should be clear that an excluded volume approach 
is by no means a complete replacement of explicit calculations of short-range correlation effects, such as through a 
Bethe-Goldstone-type of approach pp[| . There is one attempt to include short-range quark-quark correlations in the 
QMC model |26| and its further investigation in the context of a Bethe-Goldstone approach is an interesting new 
direction that should be undertaken in the near future. 
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TABLE I. Parameters used in the calculation. 



m q (MeV) 


R (fm) 


B 1/4 {MeV) 


zo 


m„(MeV) 


mu, (MeV) 





0.6 


211.3 


3.987 


550 


783 





0.8 


170.3 


3.273 


550 


783 



TABLE II. quark-sigma and omega-nucleon are used for different cases in our calculation for R — 0.6 fm. Effective nucleon 
radius, effective mass, bag eigen value and compressibility are given for different sets at the saturation density. Note that 
hard-core radius r = fm corresponds to normal QMC model 





Hard-core 


9 q s 


9oj 


R*/R 


M*/M N 


x* /x 


K 




radius (fm) 












MeV 







5.98 


8.95 


0.9934 


0.7757 


0.8659 


257 




0.4 


5.93 


8.81 


0.9936 


0.7789 


0.8684 


285 


B— const ant 


0.5 


5.87 


8.66 


0.9939 


0.7824 


0.8711 


316 




0.6 


5.76 


8.38 


0.9942 


0.7887 


0.8759 


372 







4.32 


9.87 


1.0849 


0.7388 


0.8882 


268 




0.4 


4.26 


9.72 


1.0844 


0.7426 


0.8909 


297 


B — B* 


0.5 


4.20 


9.57 


1.0839 


0.7468 


0.8938 


330 




0.6 


4.09 


9.29 


1.0830 


0.7543 


0.8988 


386 



TABLE III. Same as Table II for R = 0.8 fm 





Hard-core 


gi 


g^ 


R*/R 


M*/M N 


x* jx 


K 




radius (fm) 












MeV 







5.74 


8.19 


0.9930 


0.8034 


0.8342 


249 




0.4 


5.69 


8.06 


0.9932 


0.8060 


0.8371 


277 


B=constant 


0.5 


5.64 


7.91 


0.9935 


0.8088 


0.8404 


307 




0.6 


5.54 


7.64 


0.9938 


0.8139 


0.8461 


361 







4.14 


9.34 


1.0799 


0.7609 


0.8594 


261 




0.4 


4.09 


9.20 


1.0795 


0.7640 


0.8627 


290 


B — B* 


0.5 


4.03 


9.05 


1.0792 


0.7675 


0.8661 


322 




0.6 


3.93 


8.79 


1.0785 


0.7737 


0.8723 


378 
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FIG. 1. The energy per nucleon of nuclear matter as a function of p/po for different hard-core radii. All curves are for the 
same set of quark-meson coupling constants. The upper pannel of the figure corresponds to R — 0.8 fm and the lower one is 
for R = 0.6 fm. 
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FIG. 2. The energy per nucleon of nuclear matter as a function of p/po for different hard-core radii. The quark-meson 
coupling constants are refitted such as to obtain the correct saturation point. The upper pannel of the figure corresponds to 
R = 0.8 fm and the lower one is for R = 0.6 fm. 
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FIG. 3. The pressure of the nuclear matter as a function of the energy density corresponding to the different hard-core radii. 
The upper pannel of the figure corresponds to R — 0.8 fm and the lower one is for R — 0.6 fm. 
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FIG. 4. The effective radius of the nucleon as a function of p/po corresponding to the different hard-core radii. The upper 
pannel of the figure corresponds to R — 0.8 fm and the lower one is for R = 0.6 fm. 
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FIG. 5. The in-medium nucleon mass as a function of p/po corresponding to the different hard-core radii. The upper pannel 
of the figure corresponds to R = 0.8 fm and the lower one is for R — 0.6 fm. 
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